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The expansion of the fireball created in Au-Au collisions at ^/sjvjv = 200GeV is described in 
(3 + l)-dimensional viscous hydrodynamics with shear and bulk viscosities. We present results for 
the transverse momentum spectra, the directed and elliptic flow and the interferometry radii. 
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I. INTRODUCTION 

Heavy-ion collisions at ultrarelativistic energies per- 
formed at the BNL Relativistic Heavy Ion Collider 
(RHIC) and the CERN Large Hadron Collider (LHC) Q 
have shown that dense matter is formed in the interac- 
tion region. The fireball expands and a sizable collective 
flow develops. Effects of the flow are observed in particle 
spectra, elliptic flow, and interferometry radii. Nuclear 
modification of high p± particle spectra is understood as 
the energy loss of partons in the dense medium. 

The dynamics of the dense and hot matter can be 
quantitatively described in terms of relativistic hydro- 
dynamics Q • A further refinement of the hydrodynamic 
approach involves a finite shear viscosity of the fluid [H-Q . 
Finite shear viscosity reduces the elliptic flow in the sys- 
tem. The comparison of the experimental data to model 
predictions for the elliptic flow could be used to estimate 
the value of the shear viscosity coefficient. Most of the 
relativistic viscous hydrodynamic calculations for heavy- 
ion collisions are done in (2 + l)-dimensions [(2 + 1)-D]. 
Such a simplification requires the assumption of boost- 
invariance of the matter created in the collision. Exper- 
imental data on particle spectra at RHIC show that no 
boost-invariant region is formed, even for central rapidi- 
ties @. Only recently, first results from a full (3 + 1)-D 
viscous hydrodynamic code have become available [8|. 

We present the results of a relativistic viscous (3 + 1)- 
D code with shear and bulk viscosities applied to Au-Au 
collisions at ^/sWn — 200GeV. Hydrodynamic calcula- 
tions are preformed starting from Glauber model initial 
conditions, with the freeze-out at 135MeV and subse- 
quent resonance decays. The use of a realistic bulk vis- 
cosity in the hadronic phase allows to lower the accept- 
able freeze-out temperature, improving the agreement 
of the spectra of pions, kaons and protons and of the 
Hanbury Brown-Twiss (HBT) correlation radii with the 
data. A low value of the shear viscosity to entropy ra- 
tio r]/s = 0.08 is consistent with the observed elliptic 
flow. The expansion with finite viscosity yields HBT 
radii closer to the data than from ideal fluid hydrody- 



namics. We present first results on the directed flow in 
(3 + 1)-D viscous hydrodynamics. 



II. VISCOUS HYDRODYNAMICS 

The relativistic second order viscous hydrodynamics 
(Toj is based on the extension of the energy-momentum 
tensor of the perfect fluid 



IT = (e+pK«"-p*r 



(2.1) 



by the stress corrections from shear tt and bulk n vis- 
cosities 



(2.2) 



The fluid energy density, pressure and four velocity are 
denoted by e, p and m m respectively. The viscous correc- 
tions are solutions of the dynamical equations 

A^AWVcV^ = _H -n^d a u a (2.3) 
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and V M = A^ v d v . We take for the relaxation time tv = 
and assume rn = t„. 

Hydrodynamic simulations show that the average value 
oirj/s must be small in order to describe the experimental 
data on elliptic and triangular flows 0, H, [HI H3 • The 
extracted value is close to the conjectured lower limit 
rj/s = 0.08 [lH if the Glauber model is used for the initial 
profile of the fireball. 

The shear viscosity coefficient to entropy ratio is not 
constant in our default calculations. One expects signif- 
icant dissipation and effective viscosity in the hadronic 
cascade in the last stage of the collision 0, [l5| . In the 
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FIG. 1. (Color online) Temperature dependence of the ratio 
of shear and bulk viscosities to the entropy. The solid line 
represents the default value, corresponding to 77/s = 0.08 in 
the QGP phase and increasing in the hadronic phase to 0.26 
at the freeze-out, the dashed line corresponds to a constant 
value of rj/s — 0.08, and the dashed-dotted line represents the 
bulk viscosity f/s = 0.04 in the hadronic phase. 



hydrodynamic model, without a hadronic cascade after- 
burner, it would mean that in the hadronic phase 77/s 
increases 0, EH . The viscosity to entropy ratio is taken 
in the form 



V) 

s 



—fH G (T) + (l-f„ G (T))^^ (2.6) 



with r/ HG /s = 0.5, r) QGP /s = 0.08, and Jhg{T) = 
l/(exp((T-Tff G )/AT) + l), where T HG = 130MeV, 
AT = 30MeV. Bulk viscosity is expected to be negligible 
in the high temperature plasma phase and it must be fi- 
nite in the interacting gas of massive hadrons [T3|. We 
put a nonzero bulk viscosity coefficient in the hadronic 
phase 



(T) 
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(2.7) 



with (hg/s = 0.04 and f c {T) 
l/(exp((T-r c )/AT c ) + l), where T c = 160MeV, 
ATq — 4MeV. A similar value of the bulk viscosity and 
of its temperature dependence has been used in the 
description of the RHIC and LHC data with (2 + 1)-D 
viscous hydrodynamics @, flU fl9| . The temperature 
dependence of the viscosity coefficients is shown in Fig. 
[TJ At the freeze-out temperature Tf = 135MeV we have 
77/s = 0.26 and £/s = 0.04. 

The equation of state relating the thermodynamical 
quantities is a necessary ingredient for the hydrodynam- 
ical evolution. In the recent years, it became custom- 
ary for the hydrodynamical calculations to use for the 
equation of state a parametrization of the lattice QCD 



FIG. 2. (Color online) Temperature dependence of the veloc- 
ity of sound squared used in the hydrodynamic calculations, 
compared to lattice QCD results of [24| . 



data combined with a non interacting, hadron resonance 
gas model at lower temperatures [3, [20l - [23l |. Such an 
equation of state with a cross-over transition from the 
plasma to the hadronic phase yields a much better de- 
scription of the measured interferometry radii. The reso- 
lution of the HBT puzzle is a strong argument in favor of 
the present quantitative understanding of the dense mat- 
ter equation of state at zero baryon density [U, [22j ■ In 
the present paper we follow the prescription of Ref. [2(| , 
connecting the velocity of sound in the hadron gas below 
I45MeV to lattice QCD values above 175MeV. The in- 
terpolation between the two limiting forms is such that 
the entropy from lattice QCD is reproduced at high tem- 
peratures [2fj|. We use the recent lattice QCD results 
of the Wuppertal-Budapest group 24J. The velocity of 
sound and pressure as function of temperature are shown 
in Figs. H]and[31 

The initial density profile p(rj\\ , x, y) for the hydrody- 
namic evolution in the space-time rapidity T)\\ and the 
transverse plane (x, y) at an impact parameter b is taken 
in the form 



p(b,V\\,x,y) = 



(y b + f lll )N + + (y b -r lll )N^ 
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and the entropy density is 



i(v\\, x ,y) = «o 



p( b ,V\\,x,y) 
p(0,0,0,0) 



/fall) (2-8) 



(2.9) 



where the density in the transverse plane is proportional 



to a combination of participant nucleon p pa 



N+ + N_ 



and binary collision p bin densities, with a — 0.125. The 
densities of the right and left going participant nucleons 
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FIG. 3. (Color online) Temperature dependence of the pres- 
sure P/T 4 , compared to lattice QCD results of [23 ]. 



N±(x,y) are calculated from the Glauber model. The 
first factor on the right hand side of (I2.8I) gives a tilt of 
the source away from the beam axis (y& is the beam ra- 
pidity). It is motivated by the observed asymmetry of 
the emission of the participant nucleons. A participant 
nucleon emits particles preferentially in the same rapid- 
ity hemisphere 25 1. The hydrodynamic evolution of a 
tilted source generates the directed flow of particles, as 
observed in Au-Au collisions at ^/smm — 200GeV 



y/s^N = 200GeV [2jj. 
In the viscous hydrodynamic calculation we use the par- 
ticipant eccentricity for the initial fireball [13] ■ The op- 
tical Glauber model gives the standard eccentricity. We 
rescale the density in the transverse plane (x — > x//3, 
y — > y/3) to reproduce the participant eccentricity ob- 
tained from a Glauber Monte-Carlo model 28] , with 
/3 = 1.02-1.03. The parameters and the longitudinal pro- 
file 







(2.12) 



together with the equations for the stress corrections 
(|2.3[> and (I2.4[) are solved numerically in the x, y, m co- 
ordinates starting from tq. The initial flow is the Bjorken 
scaling flow u M = (t/r, 0, 0, z/r), the initial shear stress 
tensor takes the Navier-Stokes form, and II(to) = 0. The 
evolution in r = y/t 2 — z 2 is performed in a two-step 
method with spatial gradients calculated on a grid with 
spacing Ax = Ay = 0.24fm. At small times, viscosity 
corrections to the pressure are substantial. The formal- 
ism of second order viscous hydrodynamics is not appli- 
cable in that case [29]. To regularize we use the 
formula 



•reg 
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1/4 ■ 



(2.13) 



assuring that the longitudinal pressure does not get neg- 
ative, even in the early phase of the evolution. 

At the freeze-out temperature of 135MeV and 145MeV 
for viscous and ideal fluid hydrodynamics, particles are 
emitted from the freeze-out hypersurface according to the 
Cooper-Frye formula. Viscous corrections to the equilib- 
rium momentum distribution /o 



/ — /o + S ,/ 'shear + <5 fbu. 
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(2.14) 



yield a change in the energy-momentum tensor in 
hadronic phase 



(27T 



PV (/o + <*/) = ?T + 5T^ 



fE 1 

the sum is over all the hadron species. The corrections to 
the energy-momentum tensor fulfill the Landau matching 
conditions 



f(v\\)=^P[- 2 J fl(hl-ttO) (2-10) 

are adjusted to reproduce the charged particle distribu- 
tion in pseudorapidity. We have r/ = 1.5, er j; = 1.4 for 
the viscous evolution and 770 = 1.7, a v = 1.4 for the per- 
fect fluid evolution. The parameters of the Woods-Saxon 
density distribution in Au nuclei 

p{x,y,z) = — z r , 

1 + cxp K Vx 2 + y 2 + z 2 - R A )/aj 

(2.11) 

are po — 0.17fm , Ra = 6.38fm and a = 0.535fm; and 
the nucleon-nucleon cross section is 42mb. For the hy- 
drodynamic evolution starting at r = 0.6fm/c, the max- 
imal entropy density sq corresponds to a temperature of 
380MeV, i.e. an energy density of 33GeV/fm 3 . 



u^ST^Uu = (2.16) 

and 

u„SN£ = (2.17) 

where 

n v ' 

is the change in the conserved charge bk (e.g. baryon 
number, strangeness) in the system. 

The form of the stress corrections to the energy- 
momentum tensor and the matching conditions do not 
determine uniquely the nonequilibrium corrections 5f, 
either in the form of the momentum dependence or the 
contribution of different hadrons in a multicomponent 
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system. We use a quadratic form for the shear viscosity 
corrections 



5 f shear — /o (1 ± /o) 



1 



2T 2 {e + p) 



PW, (2.19) 



and an asymptotically linear form for the bulk viscosity, 
based on the relaxation time approximation 0, H3] , 



5 fbulk = Cbulkfo (1 ± fo) ( c^u^p^ - 



n 



with, in the local rest frame, 



(2.20) 
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d 3 p 



(2tt) 3 £ 



■/o (1 ± /o) 



s 3E 



(2.21) 

The form of the shear viscosity corrections is stan- 
dard and commonly used [31|, on the other hand, dif- 
ferent expressions for the bulk viscosity corrections are 
considered, Grad expansion |3p, e xponential (33[, and 
the relaxation time formula [7|, 130 . l34j . The assumed 
form and the relative contribution of different hadron 
species to the bulk viscosity corrections are important 
as they influence the transverse momentum spectra of 
produced particles and their relative yields. In general 
the constraints imposed by the Landau matching condi- 
tions lead to chemical nonequilibrium corrections from 
the bulk viscosity. It is true for any form of the assumed 
momentum dependence of the bulk viscosity corrections, 
linear, quadratic or exponential. Imposing chemical equi- 
librium would require Landau matching conditions not 
for the conserved quantum numbers, but separately for 
each particle specie. The actual choice of the bulk viscos- 
ity correction to be used needs a specific assumption on 
the reequilibration rates for different particles. The pro- 
posed forms of the bulk viscosity corrections range from 
assuming a similar form for all hadrons 0, HO, H^l , to a 
different form for mesons, baryons and strange particles 
34 1, to requiring the particle numbers to be unchanged 
33j. We use formula (|2.2ip assuming a common relax- 
ation time for all the particles, which leads to deviations 
from chemical equilibrium due to bulk viscosity correc- 
tions. It is a minimal assumption, but more elaborate 
ansatze are possible with different relaxation times for 
different particles for both bulk and shear viscosity cor- 
rections 34, 35]. The modification of the momentum dis- 
tribution (12.211) fulfills the Landau matching conditions 
due to the relation 
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d 3 p 



Cbulkfo (1 ± fo) 
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(2.22) 



To illustrate the effects of chemical nonequilibrium, 
let us consider a system undergoing a fast, Hubble- 
like expansion with a collective velocity u M = x M /r3, 
t 3 = y/t 2 - x 2 - y 2 



z 2 . If at some proper time T3 = to 



the interactions are turned off, the particle distributions 
at later times are 



f(p,x) = foWm 2 + ( P T 3 /t ) 2 /T dec ) 



(2.23) 



if at T3 = to the momentum distributions fo(E /Td ec ) are 
in equilibrium at the temperature Td ec - The distributions 
(|2.23|) are solutions of the Vlasov equation 



p^f(p,x) = 0. 
At times T3 > to the energy density 



(2.24) 



e(r 3 ) = E/ J^Ef n (p,T 3 ) (2.25) 



drops. In a real system rescatterings after to are still 
present and the particle distribution is driven toward the 
equilibrium with the temperature T eq corresponding to 
the energy density e eq = £(73). This equilibration is in- 
complete, and one can use an ansatz for the distribution 
with bulk corrections at the freeze-out of the form 



fo + Sfbuik = fo(s/ m 2 +p 2 X/T e ff) 



(2.26) 



with the parameters T e // and A adjusted to reproduce 
the matching conditions 



E 



d 3 p 

(27T) 5 



Ef (Vm 2 +p 2 X/T eff ) 



(2.27) 



and 



d 3 p p 2 



Ey !2;i- i - i E :-; 



f (Vm 2 +p 2 X/T eff ) = p eq + n 



(2.28) 

The freeze-out temperature T eq corresponds to the en- 
ergy density e eq , but the particle ratios correspond to the 
temperature T e ff withT e? < T e ff < T dec . The equilibra- 
tion processes drive the momentum distribution function 
from the distribution (|2.23|) towards the equilibrium one 
with the temperature T eq . The ansatz (|2.26p describes 
this effect of partial chemical and kinetic reequilibration 
after to. Using the simple two-parameter formula (|2.26j) 
the effect of the deviations from equilibrium are taken 
into account. If the chemical reequilibration processes 
are significantly slower than the kinetic ones, the particle 
ratios get fixed at some chemical freeze-out temperature 
T c h with T e ff < T c h < ?dec- Even in that case using 
T e f f , instead of forcing the particle ratios to be fixed at 
the temperature given by the energy density T eq , reduces 
the error. 

The difference between T e // and the true chemical 
freeze-out temperature T c h is not big if reequilibration 
processes are defined by the energy scales, which means 
that it is as difficult to repopulate a pion state with mo- 
mentum 800MeV as an ui state with momentum 220MeV. 
The momentum distribution with bulk viscosity correc- 
tions (|2.26p is in chemical equilibrium at the temperature 
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T e ff. But because we compare it to the reference equilib- 
rium distributions corresponding to the temperature T eq 
given by the energy density, the particle ratios are off 
equilibrium for T eqi the reason is that equilibration pro- 
cesses are not fast enough to repopulate high momentum 
states and depopulate high mass states relative to the 
instantaneous, approximate equilibrium state defined by 
the Landau matching condition. Extensive calculation in 
(2 + 1)-D viscous hydrodynamics with the bulk viscosity 
corrections of the form (|2.26l) or (|2.2ip give almost indis- 
tinguishable results for the final spectra [HI, [III [36| . In 
both cases, due to the shift in the temperature from T eq 
to T e f f the particle ratios appear as off chemical equilib- 
rium for the freeze-out temperature T eq . 

Particle emission and resonance decay is performed us- 
ing the Monte-Carlo generator THERMINATOR2 [37|. 
The hydrodynamic expansion is done using the equation 
of state at zero baryon density. At RHIC in the central 
rapidity region, the baryon chemical potential is nonzero 
yielding for the ratio of antiprotons to protons ~ 0.8. We 
reintroduce the nonzero chemical potentials at the freeze- 
out with the ratio /i/T taken from thermal model fits [38| . 
This procedure violates the baryon number flow at the 
freeze-out hypersurface, and approximately is equivalent 
to multiplying the final proton spectra by exp(/z/T) and 
the antiproton spectra by exp(— n/T). The justification 
for this procedure is that the equation of state is expected 
be moderately changing with fi at small baryon densities 
[39l | and that the energy is conserved at the freeze-out 
to the order /x 2 /T 2 . The net effect is mainly the rescal- 
ing of the relative numbers of protons and antiprotons, 
which is crucial for comparing with experimental spectra 
at central rapidities. 

III. RESULTS 

The distribution of charged particles in pseudorapidity 
is shown in Fig. 2] for different centralities. The width of 
the initial distribution of matter for the hydrodynamic 
evolution in Eq. (|2.8p is adjusted to reproduce the final 
charged hadron distribution. It is interesting to compare 
the parameters for the viscous and perfect fluid evolu- 
tions. The initial width for the viscous hydrodynamics 
is smaller. A similar behavior of the matter distribution 
in the longitudinal direction in the (3 + 1)-D evolution 
has been observed in Ref. 40]. It is contrary to the 
expectations from simple 1 + ID viscous hydrodynamic 
calculations fill ]. The reduced longitudinal pressure in 
the initial stage of the evolution 

P + 7T ZZ (3.1) 

should lead to a reduced expansion of the matter in space- 
time rapidity. Such a reduced expansion is observed in 
1 + ID calculations, using narrow initial distributions in 
space- time rapidity [Il|,[42j]. In (3 + 1)-D evolution, the 
initial distribution has a broad plateau in space-time ra- 
pidity, where no expansion occurs at early times. In the 




FIG. 4. (Color online) Pseudorapidity distribution of charged 
hadrons for centrality classes — 6%, 6 — 15%, 15 — 25%, 
25 — 35%, 35 — 45% and 45 — 55% calculated in the viscous 
and perfect fluid hydrodynamics (solid and dashed lines re- 
spectively) compared to PHOBOS Collaboration data (dots) 
[43| . The squares represent the BRAHMS Collaboration data 
for centrality — 5% 0. 

tails, outside of the plateau, the expansion is faster in the 
ideal fluid case (see the tails of the distributions in Fig. 

a. 

The centrality dependence of dN/drjps is reproduced 
using the initial entropy density scaled with a combina- 
tion of participant nucleons and binary collisions. The 
parameter a = 0.125 for the admixture of binary colli- 
sions is smaller than seen in the final density at rjps = 
(a = 0.145 [46]). The difference comes from the interplay 
of the longitudinal and transverse expansions, and from 
the entropy production in the viscous hydrodynamics. 

Pion spectra in transverse momentum are well repro- 
duced at different centralities (Fig. [S]), for p± < 1.2GeV. 
The role of the bulk viscosity is essential in reproducing 
the spectra, as it reduces the effective thermal motion of 
the emitted particles. The collective component is larger, 
corresponding to a lower freeze-out temperature. Very 
similar results are obtained using ideal hydrodynamics, 
but at a higher freeze-out temperature 145MeV. It means 
that the p± distributions in the perfect fluid case are ob- 
tained using a smaller collective flow, but larger thermal 
motion. 

Kaon spectra are well reproduced in central collisions 
(Fig. [5]). In semi-peripheral collisions the number of 
kaons is overpredicted. It may be a sign of the incom- 
plete equilibration of strangeness in peripheral collisions 
[48| . The same as for pions, the perfect and viscous fluid 
hydrodynamics give similar results. Hydrodynamic cal- 
culations describe the proton spectra up to p± < 2GeV 
(Fig. [7]) . Small differences can be observed between per- 
fect fluid and viscous calculations. Viscosity leads to 
harder spectra for protons, as heavy particles are more 
sensitive to the collective flow. Another difference is that 
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FIG. 5. (Color online) ty + transverse momentum spectra for 
centralities 0-5%, 5-10%, 10-15%, 15-20%, 20-30%, 30-40% 
and 40-50% (successively scaled down by powers of 0.1) from 
viscous hydrodynamic calculations (solid lines). The dashed 
lines for the centralities — 5% and 20 — 30% represent the 
results of the perfect fluid hydrodynamics. Data are from the 
PHENIX collaboration El. 



the number of protons is larger in the viscous calcula- 
tions, although the freeze-out temperature is lower, and 
one would expect a smaller thermal rate of production. 
This is an effect of the bulk viscosity, that drives the sys- 
tem out of chemical equilibrium. In an expanding system 
with bulk viscosity the ratio of the number of heavy to 
light particles is larger than predicted in chemical equi- 
librium at Tf (Sect. HI)) . The spectra for pions, kaons, 
and protons are very similar as obtained in (2 + 1)-D 
viscous hydrodynamics with bulk and shear viscosity 0] ■ 
A characteristic, for which viscous evolution in (3 + 1)- 
D may be important is the pseudorapidity dependence of 
the elliptic flow. Ideal fluid hydrodynamics gives a flat 
dependence, unlike measured in experiments. Dissipa- 
tivc effects in the hadronic cascade bring the results of 
the simulations closer to the data [T3|. In terms of the 
viscous hydrodynamics, one expects stronger shear vis- 
cosity corrections at forward rapidities, where the matter 
freezes-out earlier |52j. (3 + 1)-D viscous hydrodynamic 
calculations by Schenke et al. give a flat dependence of v 2 
on the pseudorapidity @, H(j , both using the average and 
fluctuating initial conditions. The same can be observed 
in our calculation (Fig. [8]). We study the effect of the in- 
crease of the shear viscosity in the hadronic phase (solid 



FIG. 6. (Color online) The same as Fig but for K + . 
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FIG. 7. (Color online) The same as Fig [S] but for protons. 
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FIG. 8. (Color online) Pseudorapidity dependence of the 
elliptic flow coefficient for charged particles for centralities 
15 — 25%, for the viscous hydrodynamic expansion with in- 
creasing (solid line) and constant (dashed line) shear viscosity 
to entropy in the hadronic phase, data from the PHOBOS col- 
laboration are denoted by dots [53]. 
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FIG. 9. (Color online) Elliptic flow of charged particles as 
function of transverse momentum from viscous hydrodynamic 
calculations for centralities c = 10 — 20% (lower line) and 
c = 30 — 40% (upper line) compared to PHENIX collaboration 
data El. 
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FIG. 10. (Color online) Elliptic flow of charged particles as 
function of transverse momentum from viscous hydrodynamic 
calculations for centralities c = 20—30% (solid line), from per- 
fect fluid hydrodynamics (dashed-dotted line), from the vis- 
cous hydrodynamics without S f 'shear and Sfbuik corrections 
at the freeze-out (dashed line), and from the viscous hydro- 
dynamics without Sfbuik corrections at the freeze-out (dotted 
line). Data are from the PHENIX collaboration [4!|. 
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FIG. 11. (Color online) Directed flow in Au-Au collisions, per- 
fect fluid (dotted line) and viscous fluid hydrodynamics (solid 
line) are compared to experimental data from the PHOBOS 
and STAR collaborations [50l. l5l|. 



line) as compared to a calculation using a constant r\j s 
(dashed line). We observe a minor improvement of the 
agreement with the data when the viscosity increases, 
as would have been expected if the effect determining 
the shape of vi(r\ps) were the dissipation in the hadron 
cascade; it is not enough to remove the discrepancy with 
the PHOBOS measurements. In the simulation with con- 
stant r\j s we start the evolution with the standard eccen- 
tricity, given by the optical Glauber model. When using 
fluctuating initial conditions, the initial eccentricity is 
the participant eccentricity, but the expansion of lumpy 
initial conditions reduces the final flow @. These effects 
depend on the centrality, the coarse-graining of the initial 



fluctuations and on the viscosity [40|, |53[ . 

The elliptic flow of charged particles as function of p± 
is shown in Fig. [9] and [10] for three different centrality 
classes. Viscous hydrodynamics gives a satisfactory de- 
scription for p± < 1.5GeV. The reduction of the elliptic 
flow from viscosity happens in the hydrodynamic phase 
(dashed-dotted line versus dashed line in Fig. [TU|) . An 
additional reduction of the azimuthal asymmetry hap- 
pens due to the shear viscosity corrections at freeze-out 
(dotted line in Fig. [TU)) . The inclusion of bulk viscosity 
corrections increases the elliptic flow slightly, as noted in 
[32| . The reason is that bulk viscosity reduces the ther- 
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FIG. 12. (Color online) HBT radii for Au-Au collisions at 
centralities — 5%. Ideal fluid results (dashed lines), viscous 
hydro-dynamic results (solid lines), and STAR Collaboration 
data [54 ] (squares) are shown. 



mal motion of the emitted particles and the momenta 
are more aligned with the collective flow of the fluid. 
The fluid velocity is transverse to the shear stress tensor 
w M 7r M „ = 0, so the shear viscosity correction in Eq. (I2.19[) 
is reduced. 

The expansion of the tilted source (Eq. 12.81) gives a siz- 
able negative directed flow (Fig. [TT]). The perfect fluid 
dynamics gives a larger v\ , as predicted in Ref. [55| . The 
formation of the directed flow from the tilted source, in- 
volves the simultaneous action of the transverse and lon- 
gitudinal pressures in the fluid (26j and it happens early 
in the evolution. Shear viscosity corrections reduce the 
longitudinal pressure (Eq. 13.11) and increase the trans- 
verse one, as a result less directed flow is generated. The 
directed flow observable is potentially a very sensitive 
measure of the pressure imbalance at the early stage. 
However, significant uncertainties are still present related 
to the initial conditions, the starting time of the evolu- 
tion, and the nature of the initial pressure imbalance [56| . 

The HBT radii arc calculated from the correlations of 



identical pions. Correlated pairs are generated from the 
Monte-Carlo events [13] ■ The interferometry radii as a 
function of the momentum pair k± are shown in Fig. [12] 
The calculations are in good agreement with STAR data; 
with viscous hydrodynamics being slightly better. The 
size of the emitting source R S ide is correctly predicted. 
The larger collective flow in the calculation involving 
both the shear and bulk viscosities yields a smaller value 
of the ratio R out /R S ide- We find the role of bulk viscos- 
ity important to obtain a common description of both the 
HBT radii and the transverse momentum spectra. The 
results for the HBT radii in the (3+ 1)-D calculations are 
very similar as in (2 + 1)-D for both the viscous @, Gil 
and perfect fluid [57| dynamics. 



IV. CONCLUSIONS 

We present a set of (3+l)-D viscous hydrodynamic cal- 
culations for Au-Au collisions at ^/snn = 200GeV. Using 
an independently developed hydrodynamic code coupled 
to a statistical emission and resonance decay event gen- 
erator [37| we evaluate several soft observables for final 
particles. The measured spectra of pions, kaons, and 
protons are well reproduced. The elliptic flow of charged 
particles as a function of p± is calculated for different cen- 
tralities. A small value of the shear viscosity rj/s — 0.08 
for calculations using Glauber model initial conditions 
yields results compatible with the data. The dependence 
of the flow coefficient vi on pseudorapidity is too flat, 
even when taking into account the increase of the shear 
viscosity in the hadronic phase. The directed flow of 
charged particles is calculated for the viscous fluid evo- 
lution, finding a reduction compared to the perfect fluid 
case. The interferometry radii as a function of the mo- 
mentum of the pair of pions are calculated. Viscosity im- 
proves the agreement with the data, especially for R S ide 
and the ratio R ou t/Rside- 

The values for some of the observables studied, the 
transverse momentum spectra, the elliptic flow as a func- 
tion of p± and the HBT radii are similar as in (2 + 1)-D 
viscous hydrodynamic simulations. The elliptic flow as a 
function of pseudorapidity is similar to previous (3 + 1)- 
D calculations using averaged initial conditions p| |4C| . 
The reduction of the directed flow from viscosity comes 
as expected [HEf . 
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